Low-Energy Structure of Heisenberg Ferrimagnetic Spin Chains 
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Static and dynamic structure factors of Heisenberg ferrimagnetic spin chains are numerically 
investigated. There exist two distinct branches of elementary excitations, which exhibit ferromag- 
netic and antiferromagnetic aspects. The ferromagnetic feature is smeared out with the increase 
of temperature, whereas the antiferromagnetic one persists up to higher temperatures. The scat- 
tering intensity is remarkably large at lower boundaries of the ferromagnetic and antiferromag- 
netic spectra. All these observations are consistent with the ferromagnetic-to-antiferromagnetic 
crossover in the thermal behavior which has recently been reported. 
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Recently considerable attention has been directed to 
quantum magnets with two kinds of antiferromagneti- 
cally exchange-coupled centers. Since coexistent spins 
of different kinds do not necessarily result in magnetic 
ground states, it remains a stimulating problem whether 
the system is massive or massless. For example, an align- 
ment of elementary cells with spins 1 , f , 1/2, and 1/2 in 
this order results in a singlet ground state and an£s 
tion gap immediately above it.EP Several authorsat 
have generally discussed what is the criterion for the 
massive phases via the nonlinear-er-model technique. On 

simple alternating-spin chains 
with two kinds of spins S and 
s which are described by the Hamiltonian 



th?. Other band t h e most 
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show ferrimagnctism instead of antifcrromagnetism. 
Considering that._all the mixed-spin-chain compounds 
synthesized so farliZP exhibit ferrimagnetic ground states, 
the study on the model ([!]), the simplest example of a 
quantum ferrimagnet, is .all the more important. Al- 
though pioneering authorsB have already discussed the- 
oretically quantum ferrimagnets in the 1980s, their inter- 
est was, for instance, the quantitative dependence of the 
thermal behavior on constituent spin quantum numbers. 
The exact diagonalization of short chains was not nec- 
essarily conclusive in understanding the universal quan- 
tum ferrimagnetic bpba.vior p.snpria.lly at low tempera- 

tures. Recent argumentflBiB§BBB0B0 have 
renewed the interest in this fascinating subject .and have 
even motivated brand-new NMR measurements!^ which 
are currently in progress. 

Let us consider the Hamiltonian (|l|) under the peri- 
odic boundary condition. We set the number of unit 
cells to ALand assume that S > s. Then the Lieb-Mattis 
theoremEJ immediately shows that the Hamiltonian (Fy) 



has ferrimagnetic ground states of spin (5 — s)N. As we 
may generally expect gapless excitations from magnetic 
ground states, we here take little interest in the naive 
problem of whether the spectrum is gapped or gapless. 
Actually, gapless excitation branchesElP lie in the sub- 
spaces whose total spin is less than (S — s)N, showing 
quadratic dispersion relations at small momenta. Thus 
the low-lying excitations of quantum ferrimagnets have_a 
ferromagnetic aspect. The other excitation branchesEJ' 
lying in the subspaces whose total spin is greater than 
(S — s)N are all. gapped and reminiscent of gapped an- 
tiferromagnetsO 

The thus-revealed low-energy structure may be recog- 
nized as a combination of ferromagnetic and antiferro- 
magnetic featuresES' .and is well explained in terms of 
the spin- wave theoryfl'LP For the sake of argument, we 
briefly review the spin- wave treatment. We introduce 
the bosonic operators for the spin deviation from a Neel 
state with M = £,-(Sf + sf) = (S - s)N as 



S+ = V2S aj , S| = S ■ 



at a j , 



(2) 



Setting the lattice constant to a, we define the momen- 
tum representation of the bosonic operators as 

1 _Y^2iafc(j-l/4) 

J, j , 

(3) 
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The Bogohubov transformation 



N 

3 



-2iafc0 + l/4) & . 



a k = coshf^. a k + smh8 k b\ , 



— sinh#ju al + cosh6*t. b 



k "k 



with 



tanh20jfc = — — cos(afc) , 



(4) 



(5) 
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results in the two distinct dispersion relations 



generally defined as 



0J2 



u>kT (S - s) , 



where 



(S-s) 2 + 4Sssin 2 (afc) 
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The a- and /3-bosons, respectively, have the effects of 
reducing and enhancing the ground-state magnetization 
and may therefore be regarded as ferromagnetic and 
antiferromagnetic. The obtained dispersions u> k and 
u> k are indeed reminiscent of those of ferromagnets and 
gapped antiferromagnets, respectively, though the an- 
tiferromagnetic gap within the lowest-order spin-wave 
theory, uj~£=(nr-i J: is much smaller than the true value 
A = 1.759JBe) 

Two distinct features are most prominently exhibited 
in the temperature dependence of the specific heat.tiP 
At low temperatures a ferromagnetic T 1 / 2 initial behav- 
ior appears, whereas at intermediate temperatures we 
encounter a Schottky-like peak which can simply be de- 
scribed by the antiferromagnetic gap A. Motivated by 
these observations, we here investigate the low-energy 
structure of quantum ferrimagnets in detail. Although 
the low-lying antiferromagnetic excitations lie in the 
background of the ferromagnetic spectra, why does such 
a well-pronounced Schottky-like peak appear? This is 
the subject in thcjpiiiesent article. Since previous numer- 
ical investigationsBLP suggest that quantum behavior of 
the model (|l|) is qualitatively the same regardless of the 
values of S and s as long as they differ from each other, 
we restrict our argument to the case of (S, s) = (1, 1/2). 

In order to reveal how the ferromagnetic and antiferro- 
magnetic aspects are exhibited, let us investigate static 
structure factors as functions of temperature. Generally, 
the static structure factors of ferromagnetic mono-spin 
chains have their peaks at the zone center k — 0, while 
those of antiferromagnetic mono-spin chains at the zone 
boundary k — n/a. However, in the present case, the sys- 
tem is composed of two interacting sublattices and both 
ferromagnetic and antiferromagnetic peaks appear at the 
center of the reduced Brillouin zone —ixjla < k < it /2a. 
Therefore the position of the peak can not be used to 
determin its nature. In an attempt to find relevant oper- 
ators to detect the ferromagnetic and antiferromagnetic 
modes, we here rely on the spin-wave theory, namely, 
we regard the a-boson and /3-boson correlations as ferro- 
magnetic and antiferromagnetic, respectively. Keeping 
in mind that the present model is isotropic, thermal av- 
erages (a[a k ) and (filfa) can be calculated in terms of 
spin operators as 



(a\a k ) 
Wk) 



s-(k) = (oz k o-), 

S+(k) = (0+_ k O+), 



(8) 



with 



(9) 



■ v^sinhflfc s z k , 
V^coshflfc si , 

where the Fourier transforms of the spin operators are 



O k = cosh0 fe S% 
0+ = sinh0 fc S* 



S X = _L y c 2iafeO-l/4) 5 A 



1 



y^ e 2iafc(j + l/4) s A 
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and 6 k is defined in eq. (||). 

Employing a quantum Monte Carlo methodjjtij we 
have calculated the thus-defined static structure factors 
at various temperatures, which are shown in Fig. 1. 
Both structure factors S T (k) have their peaks at k = 0, 
which reflects a combination of ferromagnetic and an- 
tiferromagnetic features in this system. However, they 
show quantitatively different temperature dependences. 
While the ferromagnetic peak is rapidly smeared out as 
the temperature increases, the antiferromagnetic one still 
persists at rather high temperatures. This observation 
implies that the antiferromagnetic correlation is domi- 
nant at intermediate temperatures and thus well explains 
the well-pronounced Schottky-like peak of the specific 
heatlf 

The existence of a gapped antiferromagnetic mode 
does not necessarily result in a Schottky-like behavior of 
the specific heat, which is basically inherent in two-level 
systems. As long as the specific heat shows a prominent 
peak reflecting the gap A, the lowest- lying antiferromag- 
netic mode should be fully stressed in the spectrum. In 
this context we inquire further into dynamic structure 
factors. We introduce the ferromagnetic and antiferro- 
magnetic dynamic structure factors at T = as 

S-+{k,u) = EJ(n)^ + s-\Q)\H(u - (E n - E )) , 



respectively, where S k = S k ± iS%, s k — s k ± is 



(11) 



'fc> VI 

denotes an eigenstate of the Hamltonian (|1|) with en- 
ergy E n , and E is set to the ground-state energy. The 
ground states of quantum ferrimagnets are macroscopi- 
cally degenerate and thus the choice for |0) is not unique. 
We here assume the highest-magnetization ground state, 
namely, the lowest-energy state with M = (S — s)N, 
as the initial state |0) so that S ^(k,uj) and S A (k,ui) 
can detect the ferromagnetic and antiferromagnetic el- 
ementary excitations, respectively. S c " 7 {k,uj) (a = ±, 
a = —a) is generally expressed in terms of its corre- 
sponding Green's function as 



S aS {k, oj) = --ImG ff5 (fe,w). 
G" T<T (fc,cj) is given as a continued fraction!^ 

(ois^io) 



G CT<T (fc,u;) 



Id — Clo 
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where the coefficients a n and b n are obtained by the 
Lanczos method through a recursive relation 

n - n = U\f n ) - a n \f n ) - b 2 n \f n 
|/ )=S£|0>, b = 0, 



i). 



(14) 



with a set of orthogonal states {]/«)}■ This method 
works fairly welK'E!!' for the study of low- lying states 



r'errimagnetic spin i^nains 



of various one-dimensional magnets. 

The thus-calculated dynamic structure factors are 
shown in Fig. 2. We find qualitatively the same be- 
havior at all the chain lengths we treated and therefore 
the present observations may be regarded as the long- 
chain behavior. Tjhjs is due to the considerably small 
correlation lengthErB 1 of the system. In effect, the anti- 
ferromagnetic gap A is 1.7594J, 1.7592J, and 1.7591J 
atAr = 8,iV = 10, and TV = 12, respectively, the val- 
ues of which are almost seen to converge. The scattering 
intensity is remarkably large at the lower boundaries of 
the ferromagnetic and antiferromagnetic spectra at each 
k, as was expected. In particular the lowest-lying anti- 
ferromagnetic poles hold most of the antiferromagnetic 
scattering weight. They are never smeared out in the 
ferromagnetic spectra but are prominently present. As 
far as antiferromagnetic eigenvalues are concerned, we 
can in principle detect all of them in the subspace of 
M = (JS — s)N + 1, and therefore, their lower boundary 
seems to form a split band separated from an upper con- 
tinuum. Thus the low-energy structure is essentially de- 
scribed by these well-pronounced ferromagnetic and an- 
tiferromagnetic bands, which are almost parallel to each 
other and are separated by the gap A — 1.759J. This 
explains the reason why the specific heat is well fitted to 
the simple Schottky formO 
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at intermediate and even higher temperatures. 

In conclusion, we have numerically studied the low- 
energy structure of Heisenberg ferrimagnetic spin chains 
in an attempt to understand their novel thermal be- 
haviorE-P The two distinct static structure factors illus- 
trate why the clear ferromagnetic-to-antiferromagnetic 
crossover is observed in the temperature dependence 
of the thermodynamic quantities. As the tempera- 
ture increases, the ferromagnetic nature rapidly fades 
out, whereas the antiferromagnetic one persists up to 
higher temperatures. That is why the specific heat 
exhibits a well-pronounced Schottky-like peak and the 
susceptibility-temperature product shows an increas- 
ing behavior with temperature, in spite of the low- 
temperature ferromagnetic behavior. On the other hand, 
the dynamic structure factors clearly indicate that the 
present model effectively behaves as a two-level system 
unless the temperature is sufficiently low. The findings 
are a suggestive guide for neutrcprscattering measure- 
ments. Recently Hagiwara et alx3 have succeeded in 
obtaining the fundamental parameters of a bimetallic 
chain compound NiCu(pba)(D20)3-2D20, where pba = 
l,3-propylenebis(oxamato), which is a prototype of the 
(S,s) = (1,1/2) ferrimagnet. We hope that our calcu- 
lations will motivate further experimental study on this 
topic. 
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Fig. 1. Ferromagnetic (a) and antiferromagnetic (b) static struc- 
ture factors as functions of temperature for the chain of N = 32. 
The temperature is indicated in units of J/fee m the inset. 



Fig. 2. Ferromagnetic (x) and antiferromagnetic (+) dynamic 
structure factors at the absolute zero temperature. The lowest- 
energy five poles are shown for each of them. The intensity of 
each pole is proportional to the area of the circle. 
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Fig. 2(a) S. Yamamoto and T. Sakai 




Fig. 2(b) S. Yamamoto and T. Sakai 
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Fig. 2(c) S. Yamamoto and T. Sakai 



